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Notices of the Aeronautical Society of Great Britain. 


Members’ Subscriptions. 


Those Members who have not yet paid their subscriptions for the current year 
are reminded that these fall due on the 1st January each year. 


Wilbur Wright Lecture, 1916. 


The Fourth Wilbur Wright Memorial Lecture will be delivered on Tuesday, 
June 6th next, at three p.m., in the Theatre of the Royal Society of Arts, John 
Street, Adelphi, W.C. 

The title is ‘‘ The Life and Work of Wilbur Wright,’’ and the lecturer, Mr. 
Griffith Brewer, A.F.Aé.S. 


The Right Hon. the Lord Montagu of Beaulieu will preside. 


Tickets will shortly be available for distribution from the Offices of the 
Society, 11, Adam Street, Strand, W.C., to which address applications should 
be made. 


Edward Busk Studentship in Aeronautics. 


Applications will shortly be invited for the above Aeronautical Research 
Studentship. The approved scheme provides that the award and control of the 
Studentship are to be in the hands of a Committee of whom three members are 
appointed biennially by the Aéronautical Society, and three by King’s College, 
Cambridge (Edward Busk’s old College), with two of the Trustees of the Fund. 
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The remuneration of the Student will be at‘a rate not exceeding 4150 per annum 
for a period of one year, with an extension in special circumstances to two years. 
The Student will be elected by the Committee without examination and shall not 
have completed his 25th year on the 1st day of October preceding the election. 


The Student will be required to make reports to the Committee on his work 
from time to time, and at the close to present a full report to the Aéronautical 
Society, the University of Cambridge and the Trustees, and he may be required 
to deliver a lecture or lectures embodying the substance of the Report. 


The intention of the founders of the Studentship is that it should be primarily 
for the purpose of encouraging and helping young aeronautical engineers to make 
a good start on their career. 


Inquiries should be addressed to the Secretary, Aéronautical Society, 11, 
Adam Street, Strand, W.C., from which address application forms will shortly 
be procurable. 


Edward Busk Studentship Committee. 


Two of the Society’s three representatives on the above Committee have been 
appointed, viz. :—Mr. Horace Darwin, F.R.S., F.Aé.S. and Dr. R. T. Glazebrook, 
€.B., F.R.S., F.Aé.S. The name of the third representative will shortly be 
announced. 


Annual General Meeting, 1916. 


The Annual General Meeting of the Society will take place, as last year, 
immediately preceding the Wilbur Wright Lecture. 


AGENDA. 


To receive and approve the Report of the Council on the state of the Society, 
and the Balance Sheet of Aérial Science, Limited. 


To discuss and determine such questions as may be proposed by the voters 
relating to the affairs of the Society, and to fill the vacancies of the 
Council for the ensuing year. Any voter desirous of proposing any 
subject for discussion at the Annual General Meeting shall give notice 
in writing to the Secretary, which shall be received by him by noon on 
23rd May, 1916. 


The retiring Members of Council are :— 


Harris Bootnu, Squadron Commander F. K. McCLean, 
J. H. LEDEBOER, Major-Gen. R. M. Ruck, C.B., R.E., 
Lieut. A. R. Low, R.N.V.R., Lieut.-Col. F. H. Sykes, C.M.G., 

Dr. T. E. Stanton, F.R.S., Dr. R. MULLINEUX WALMSLEY, 


who are eligible for re-election. 


Nominations of candidates for election to the Council shall be signed by the 
voters proposing them (two voters and no more), and must be received by the 
Secretary by noon on 16th May, 1916, with an intimation in writing by the 
voters nominated that they are willing to serve. 


The Council will also ask for approval of their action in postponing, on 
account of the war, the date of the Annual General Meeting. 
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Elections. 


The following elections to the Society, some of which were, by oversight, not 
previously published, are notified :— 


Members : J. C. G. SykEs, C.B., dating from October, 1915. 
ALEXANDER WADDELL, ¥ October, 1915. 
J. A. WHITEHEAD, m October, 1915. 
Lieut. JoHN SowrEy, R.F.C. ,, January, 1916. 


Associate Members : FrRANcIS WarRD, dating from October, 1915. 
JosEPH HOOoLey, January, 1916. 
H. SKELTON, * March, 1916. 


Students: L. A. C. pE La GarpbE, dating from November, 1915, 
S. L. GREEN, . February, 1916. 
re February, 1916. 


P. €. THORNTON, 
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FOURIER’S THEOREM IN AERODYNAMICS. a 
BY PROF. G. H. BRYAN, 
WITH TABLES DRAWN UP BY CARADOG WILLIAMS, B.SC. 


PART I.—PLANE SURFACES. 


1. In ‘‘ Practical Aerodynamics,’’ Major B. Baden-Powell quotes the following 
formutez for the pressure of the air on a moving plane when expressed as a 

function of the angle of attack a, taking the pressure as unity when a = go”. 


I. Duchemin 
2sina 


I + sin?a 


II. Rayleigh (hydrodynamical) 
7 sina 


4+ 7Ssina 
IIa. Gerlach 


III. Renard 


IV. Soreau 


/ I \ 
sina |i + ————— 
I + tan? J 
V. Eiffel, 1907, 
fa/30 from a = 0 to a = 30° 
\1 from a = 30° to a = go° 


VI. Another formula frequently quoted is 


2 sin a from a = 0 to a = 30° 


I from a = 30° to a = 90° 


and a more complicated formula, depending on the aspect ratio of the plane, has 
been proposed by Soreau. 


The use of Fourier’s series affords a convenient means for standardising and 
comparing the results both of these formule and of experimental investigations. 
It might naturally have been supposed that this method would be the first to 
suggest itself to the earlier writers on the subject, but the fact that results are 
not given in the form of such a series in the books which are most read indicates 
that this method has been wholly or partially overlooked, or at any rate has not 
received so much attention as it appears to deserve. 


The object of the present paper is to illustrate the application of Fourier’s 
theorem to the various aerodynamical magnitudes which are functions of the 
angle of attack of a moving plane, viz., resultant thrust, lift, drag, co-ordinate 


VIIM 


| 
(4 + =) sina 
4+ 7sina 
sin a [2 (a 1) sin? a] 
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of centre of pressure and moment of resultant thrust about the centre of area. 

The formule quoted above afford a sufficient variety of examples of the application 

of the method to empirical formule, while the tables contained in Eiffel’s 

d ‘* Resistance of the Air and Aviation,’’ translated by Hunsaker, give a sufficient 
diversity of experimental data to exemplify the methods of determining the 
Fourier co-efficients from tabulated statistics. To calculate the expansions for 
every series that has been proposed or for every collection of experimental results 

_ that has been tabulated would of course be out of the question, and it is therefore 
intended that the data here given should be regarded as illustrating the methods 
and the general character of the results rather than as leading to anything final. 
As fresh experiments are made and fresh formule proposed, it should be easy 
to represent them in the same form. 


The advantages of the Fourier analysis are 

(1) The standardised form in which different results can be expressed, 
the only differences being in the values of the numerical coefficients. 

(2) The elasticity of the formule which can be adapted to functions of 
almost any kind whatever. : 

(3) The possibility of obtaining approximate formule by retaining only 
a few terms of the series, or of obtaining closer approximations in some 
cases by the addition of further terms. 


2. Duchemin’s formula can be most easily expanded by means of De Moivre’s 
theorem. Putting x = cosa + isina we find 


4 p q 


q—p \ x?—p x? —q 
where 
P= 3—2,/2,49 =3+ = 4,/2, Qp=l. 
Expanding the fractions in ascending powers of p or 1/q we obtain 


./2 1} + pcos 2a + cos 4a + cos 6a + -} 


Multiplying by 2sina and using the usual trigonometric formula we 
obtain 
2sina 
= (4— 2, /2) (sina + psin 3a + sin 5a, etc.) . (1) 
I + sin?a 
where 
Pp = 3—2,/2 = 0.1716 
This series, it will be seen, converges very rapidly and is therefore convenient 
for calculation. 


3- Sorean’s formula as quoted above is still simpler as it gives 
|\ 


= 1.25 sina + 0.25 . 2 
3 
1+tan“a 


sina |1 + 


Both the above formule hold good for all values of the angle a. If a is 
negative and between o and — 180° the pressure is on the opposite face of the 
plane and therefore negative. Also both expressions become equal to unity 


} : 

f 
4 
I 4 
—4 — 4x? a 
= I + (x — = q 
I + sin?a —6+%x74+ x7? xt—6x?4 1 

4x? 
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when a= 90° and therefore they correctly represent the pressure ratio 
P (a): P (90°). 

Rayleigh’s formula II. can only hold good between o and 180° and it is 
obvious from symmetry that Renard’s formula III. cannot hold good beyond 90°, 
as it would make P (180° —a) unequal to P (a), where P is the pressure. 


Rayleigh’s expression does not become equal to unity when a = go°, but 
Gerlach’s modification (III.a) does so. In the case of Renard’s formula we 
must divide by 3—4z to obtain the pressure ratio P (a): P (go°). 


In all the formule, except I. and IV., the expansion must be performed by 
the usual Fourier methods. 


4. Digression on Fourier Series.—The correct form of Fourier’s theorem 
states that if f(a) is a function of the angle a which is given for all values of a 
or} 0 to 360°, and which outside these limits satisfies the condition f (n.360° + a) 
= f(a) then f(a) can be expanded in a series of sines and cosines of multiples 
of a, and the coefficients of cos na and sin na can be obtained by multiplying by 
cos ma and sinna respectively and integrating between the limits o and 360° 
(or o and 27 since the integration usually has to be effected in radians). 


If the value of f(a) is only given for angles in the first quadrant it would 
be possible to find any number of different series representing the function, by 
assigning arbitrary values to f(a) in the remaining quadrants, but the following 
forms are the most convenient. 


A f(a) = f (180° — a) = f (180° + a) = f(—a) 


Even cosine series 
a, + a, COS 2a + a, cos 4a + 


Odd cosine series 
ad, COS a + d, COS 3a + COS 5a + 


& f (a) = — f (180° — a) = + f (180° + a) - Te 


Even sine series 
b, sin 2a + b, sin 4a + b, sin 6a + 


D f (a) = f (180° — a) = — f (180° = i a) = —f(—a) 
Odd sine series 
b, sina + b, sin 3a + b, sin 5a + 


Having regard to the conditions of symmetry, we see that if the formule 
are to hold good for angles of any magnitude in the case of a plane symmetrical 
lamina such as a rectangle— 

1. The ratio P (a): P (g0°) must be expressed in the form of a series of sines 
of odd multiples of a, or as we say for shortness, an odd sine series, and it will 
be seen that this form is uniform with the expansions of Duchemin’s and Soreau’s 
formule given above. We shall call this expansion 

P, sina + P,sin3a + P;sin5a+.... (3) 


2. The distance of the centre of pressure from the centre of figure must be 
expressed as an odd cosine series, say, 


C COSa + cos 3a + C, ,;COSSa+... (4) 


3. The moment of the resultant pressure about the centre of figure must be 


expressed as an even sine series, say, 


G, sin 2a + G, sin 4a + G, sin 6a. ‘ ‘ ‘ (5) 


XUM 
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5. Modification for Small Angles—Cosine Series.—If we use the formula 
f(a) = P, sina + P, sin 3a + P; sin 5a + 
for the ratio f (a) = P (a): P (g0°) it will be seen that when a is small the ratio 
of f(a) to sina approximates to the value 
P, + 3P, + 5P; + 


Although this series is really convergent in the cases which occur in 
aeronautics it js obviously very inconvenient for calculation as it may converge 
very slowly. To meet this difficulty it is sometimes better to expand f (a)/sina 
in the form of an even cosine series so that, say, 

f(«) = sina x (Q, + Q, cos 2a + Q,cosga + .. .) 
= Q, sina + 4Q, (sin 3a — sina) +4Q, (sin 5a — sin 3a) + 


Thus 


P, Q.— 42:, P, + (Q.— P, 4 (%,— Q,). (6) 


Knowing the cosine series for f(a)/sina the corresponding sine series for 
f(a) is thus easily deduced, but the converse is not possible unless the value of 
is known. Of course we have, however, 

Q =P, +P, +P, + 
which may be used if the series is sufficiently convergent or if a small number of 
terms give the necessary degree of approximation. 


0 


If, however, Q, has been calculated independently the remaining coefficients 
of the @ series can easily be deduced from the P series. 


6. Expansion of Eiffel’s Form (V.)—By multiplying and integrating between 
o and $a we get, using circular measure and writing x for a, 


6 | 
—P,= —xsinnxdx + | sinnx dx 
4 
Jo 7/6 
24 
P, = ——sin— 
6 
Putting n = 1, 3, 5, 7, etc., the expansion gives 
12 /sina 2 Sin 3a sin 5a sin 7a 2 sin 9a sin Ila 
f(a) = — + — + 
(7) 
the coefficients containing factors which recur in the order 1, 2, 1, —1, —2, —T1, 


1, 2, 1. It will be seen that the coefficients change sign after the first three, 
which is not the case in Duchemin’s formula. 


For the corresponding form (VI.) we find 


7 sin (n — 1) 7/6 sin (n + 1) 7/6 
—P, = + 
4 n (n — 1) n(n + 1) 
P,=—+ 
4 


| 
| 
| q 
a 
| 
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and the series becomes 


2 sina 2,/3 (sina 2 SiN 3a sin 5a 


+ + —— + 
3 7 1.2 2.4 4-5 
sin 7a 2 sin ga sin Ila sin 13a 
7.8 8.10 10.11 13.14 
while if it be required to use the cosine series for f(a)/sina we have 
2 2 log,, (2 + ./3) 
Q, 
3 
7. Renard’s formula, III., when expanded in the same way leads to 
(11 \ I 
sin a 2— (a— 1) sin?a} = — — — — | sina— — —|] sin 3a 
4 16 a] \ 4 16 
I 3 I (3 1 \ {4 I 
+—— + —| singa— |— + —]sins5a + |— + —] sin7a 
\2? \ 6? 2* 6? 107 
—|— + —]singa+ |— + sintta, ete. (9) 
\10? 6? / \10? 147] 
while the constant term of the cosine series for 12 — (a —1)sin?a } is equal to 
5 7 I 
2 8 27 


8. Rayleigh’s and Gerlach’s formule.—The determination of the coefficients 
in the expansion in this case depends on the evaluations of integrals of the 
form 

m/2 
cos nx dx sin nx dx 4 
——— or —— where p = — 
sinx + k sinx +k T 


which we will denote by J, and J,, the choice depending on whether it is decided 
to use the sine or cosine series. The evaluation of I,, J, and J, is discussed in 
all elementary text books (for example, Edwards’s Integral Calculus for 
Beginners, §§74, 75, gives all that is necessary). For the subsequent terms of 
the series we may use a formula of reduction based on the formule 
cos nx + cos (n — 4) x = 2 cos (n — 2) x (1 — 2 sin? x) 
sin nx + sin (nm — 4) x = 2 sin (n — 2) x (1 — 2 sin? x) 
Now 1 —2sin? x = 1 — 2k*— 2(sinx + k) (sinx —k) and we easily deduce 
m/2 
I, + (4k? — 2) I,_2+ In_g= —4 (sin x — k) cos (n — 2) xdx . (10) 
Jo 
and a similar form for the J, integrals. The only objection to this method is 
that since 4k? — 2 is roughly equal to 4.4, a small error in the value of one term 
will produce increasing errors in the subsequent terms, and if the calculations are 
carried far enough the effects of such an error will be to make the coefficients 
ultimately divergent as Mr. Creak has verified by actual calculation. For the 
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calculations required in aeronautical applications all that is necessary is to calculate 
the first terms to a few extra decimal places in order to get the subsequent ones 
correct. In some of Mr. Creak’s calculations an alternative formula, obtained 
by the use of De Moivre’s Theorem, has been used, and is free from this 
objection. 


The numerical coefficients in the various expansions have been evaluated by 
Mr. T. G. Creak, M.A., of Llanberis, and are given in Table A. While the 
results represent the expansions of the functions in question, it is quite evident 
that for aeronautical purposes a few figures only should be retained. 


g. Calculations from experimental data.—The methods of evaluating the 
coefficients in a Fourier expansion from statistical data have been fairly well 
investigated by mathematicians and others. In most cases it is assumed that 
the data have been tabulated at equidistant intervals, say, at intervals of S° 
or 10° and in such cases the methods are well known. But in aeronautical 
experiments, such as those of Eiffel, the intervals are different in different parts 
of the scale; for example, in one of Eiffel’s experiments with a square plate, 
the data are tabulated for angles of 5, 10, 20, 30, 35, 38, 40, 45, 60, 75, 90 
degrees. Moreover, it is very important that the method adopted should be as 
simple as possible, because if it is to be of any use it must be capable of being 
employed by practical experts who have no time for making elaborate calcula- 
tions. Also it would be quite unnecessary and indeed unadvisable to carry the 
figures to a high degree of approximation, and, indeed, it will certainly be 
necessary to cut down the superfluous places of decimals in the calculations 
tabulated in illustration of this paper. 


It thus becomes necessary to resort to approximate methods of integration, 
and after examining a number of methods I have come to the conclusion that a 
generalisation of the trapezoidal rule is the best method to adopt. In other 
words, when the values of f(a) have been plotted in a graph, the quantity to 
be expanded is taken to be represented by the polygonal outline formed by 
joining the points up by straight lines instead of a continuous curve passing 
through the points. The formula required is obtained as follows :— 


Taking 
y = f(a) = P, sina + P, sin3a + P, sin 5 
‘we have on multiplying and integrating, and then integrating by parts and 
1eplacing a by x,* 


T f (x) cos nx | f! (x) cos nx 
‘ 


dx (11) 


n 


[ 

P= f (x) sinnx dx = — | 

| lL on jo | 

Jo Jo 

In the case of pressure on a plane lamina f (0) = 0, also cosnz/2 = 0 since n 
is odd, the term in square brackets therefore vanishes. 


For the polygonal outline if (x,, y,) and (x, 41) 1) are two consecutive 
vertices we have for this element 
Yr+1 —Yr Ay 
fi (x)= = |— 


* x is taken to denote ‘the independent variable in the integration formule while a 
«denotes the angle of which the function. is expanded. 


| 
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say, and is constant in the interval. This gives 


[ Ay sin nx, 4 sin nx, 
as 
4 | Ax J, n? 
49 r Ay 
= — nx, | — 
n? \| ax Ax 
which we may write 
I Ay 
=— — XsinnxA 
n? Ax 


In this formula angles are measured in radians. In practical work it is 
necessary to measure Ax in degrees, and when this is done we find 
4 x 180 [ Ay ]r Ay ]rt+1 
P, = ———— sin nx, + | 
mn? (| Ax Ax |. 


(12) 


The values of the logarithm of 4 x 180/7*n? for different values of n have 
been tabulated in Table C. 


To illustrate the method of calculation a specimen calculation is shown in 
Table D. 


As some of the coefficients are negative and their logarithms are imaginary, 
I have used the prefix j to denote log (— 1) must be added. If 7 occurs twice in 
the logarithms to be added together the two j’s of course cancel; if j occurs in 
the final result the quantity of which it is the logarithm must have the minus 
sign prefixed. 


10. Mr. Caradog Williams has applied the method to the statistics of Eiffel’s 
experiments tabulated on pages 118, 119, of the English translation of his book, 
and the results of Mr. Williams’s calculations are tabulated in Table B. 


It will be seen that in some cases the terms after decreasing again increase, 
so that after the first few terms further convergence does not take place for some 
time. An examination of the calculations has shown us that these large terms 
are due to the discontinuity which in the square plate occurs between 35° and 40°. 
In fact in the expression for P, given by (12) the value of A (Ay/Ax) is 
large in the neighbourhood of the discontinuity and small elsewhere. Conse- 
quently the terms contributed by this part of the scale considerably exceed the 
algebraic sum of the positive and negative ones contributed by the other parts of 
the scale. 


It will be seen that the series derived from these experimental data usually 
converge more slowly than those proposed from theoretical considerations. On 
the other hand the results differ so widely in different cases that it hardly appears 
likely that they could adequately be represented by any better method than the 
Fourier expansion. In order to estimate the number of terms required to give a 
fair approximation to the correct result Mr. Caradog Williams has drawn up 
Figs. 1, 2 and 3, in which the original graph of Eiffel’s tabulated results is shown 
side by side with the curves obtained by taking a limited number of terms of 
the Fourier expansion. It will be seen that in the Fourier curves the fluctuations 
are in general rather less marked than in the original curves, the former curves 
lving as a rule on the concave side of the latter. This difference is no doubt 
in some measure due to the substitution of a polygonal outline for a continuous 
curve in the evaluation of the coefficients. This polygonal figure lies on the 
concave side of the actual curve and when the higher terms of the expansion 
are neglected, we naturally expect to get a curve lying within the original one. 
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It is obvious that other methods .of evaluating the coefficients are possible in 
which this defect is not liable to occur, but as previously stated, the calculations 
in such cases would probably be more difficult and it would hardly be worth 
while to make the attempt at present. 


On the whole I should think that in most applications it would be of little 
use to carry the calculations beyond two places of decimals and that a formula of 
the form P, sina + P, sin3a + P,sinsa will be sufficient for most practical 
purposes except in the case of the square plate and a few other instances, in 
which it may be necessary for greater accuracy to retain sin 7a, sin ga, or even 
further terms. 


It will be seen that all the formule expanded in Table A are represented 
with sufficient accuracy by the first three terms, the neglected coefficients in no 
case exceeding about 2 per cent. of the coefficient of sin a. 


11. A possible criticism of this conclusion may be anticipated. It will, no 
doubt, be said that when all is said and done we have only got a result which 
can be reduced to the form A sina + Bsin'a + Csin'a. This conclusion, how- 
ever, applies to nearly every practical application of Fourier’s series and still 
more so to other expansions used by mathematical physicists, in particular, to 
the use of spherical harmonies. It is certain that if equally good results could 
be obtained more easily by other methods we should do well to make a bonfire 
of most of the existing literature of the Fourier analysis. The important point 
is that the best results are obtained when the Fourier mill is used to grind out 
the successive terms of the expansion; when this has been done, it is immaterial 
in what form they are exhibited. The justification for this statement lies in the 
fact that Fourier’s method gives identically the same results that would be 
obtainable by the method of least squares. 

12. The expressions for the lift and drag, expressed as fractions of P (go°), 
reduce to the forms 

Lift = 4(P, + P,) sin 2a + 3(P, + P,) sin 4a + 4 (P; + P,) sin 6a+ 

Drag = 4P, + 4(P,;— P,) cos 2a + 4(P,— P,) cos 4a + 
(13) 

The lift and drag coefficients are calculated in the case of the square plate, 
the rectangle of aspect ratio 9, and rectangle of aspect ratio 1/3, and the results 
are tabulated in Table F. 


13. The constant term Q, of the cosine series for f(a)/sina does not lend 
itself to calculation by this method, for if we assume f (a) to be a linear function 
of a between two consecutive plotted points on its graph, we obtain an integral 
which cannot be evaluated in finite form. The only way, therefore, is to assume 
a polygonal outline for the graph of f(a)/sina in which case Q, is given by 
the trapezoidal rule. As, however, this is different to assuming a polygonal 
outline for the graph of f(o) it is obvious that the difference might have 1 
considerable effect on the subsequent calculation of the terms of the Q series 
from those of the P series. Mr. Caradog Williams has therefore calculated the 
cosine series in the case of the square plate and the rectangle aspect ratio 1.5, 
and it will be seen that the differences in this case are unimportant. The results 
are given in Table E. This objection does not of course apply in the case of 
the theoretical expansions of Table A. 


14. Formule for the pressure in oblique positions.—If a rectangle whose 
aspect ratio is 6 be rotated in its own plane through 90°, it will be described in 
the new position as having an aspect ratio of $. If, however, it be turned through 
an intermediate angle f, it will be in an oblique position relative to the wind, 
as exemplified by the case in which an aeroplane is struck by a side wind or is 
moving sideways or turning about a point in front of or behind its main planes. 
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From the absence of experimental data in the literature which I have consulted, 
it would appear that little has been done in investigating the pressures in these 
oblique positions. But considerations of symmetry show that the proper form of 
the Fourier series to use in such cases is 


a, + a,cos 2B + a,cos4B... 


where 8 may be called the azimuth of the wind relative to the plane. When the 
pressure function is expanded in sines of multiples of a each of the coefficients 
must be expanded in this form. The result will be a double Fourier series giving 
the pressure on the plane as a fraction of P (90°) for any direction whatever of 
the relative wind. 


For small angles of attack in which the value of P (a) is proportional to sina 
the proper form to assume will still be 


(a, + a,cos 2B + a,cos48 + .. .)sina 


If it should be possible to neglect the terms involving cos 48 and higher 
multiples, the expression for any such coefficient will assume the form 


c,cos*B + c, sin?B ; ‘ (14) 


where ¢, and c, represent the values of the corresponding coefficient in pterygoidal. 
and apheroidal aspects for the same plane. For example, if the dimensions of 
a rectangular plane are 15cm. and gocm., we should multiply the coefficients in 
the fifth row by cos?@ and those in the last row by sin*8 and add, thus giving 


(1.133 cos?B + 1.175 sin?) sina + (0.267 cos?B + 0.144 sin?) sin 3a + 


There is, however, no evidence that the coefficients above cos 28 can be 
neglected, and, on the contrary, such an assumption, when applied to a square 
lamina, would lead to the result that the pressure was the same, for a given 
angle of attack, whether the square was placed with a diagonal or with one of 
its edges perpendicular to the wind. These are questions which can only be 
determined by experiment. I should expect to find that the form (14) would 
be fairly accurate if applied to an elliptic lamina, but for a rectangle it would 
not be surprising to find discontinuities for certain values of the angle 8. For 
a square plate the series must necessarily be taken of the form 


a, + a, sin 4B + a, sin 88 + 


Further than this it is impossible to proceed without actual experimental 
evidence. 


15. Applications of the Method of Least Squares.—The well known connec- 
tion between Fourier’s Theorem and the method of Least Squares may be used 
with advantage in the approximations for aerodynamical coefficients. Taking, 
for example, an expansion 

f(a) = a, sina + a,singja + ... + a,sinna + 
the mean square of the error resulting from taking the first n terms of the series 
as an approximation for f(a) in the interval from o to 42 is 
2 (x/2 
(f (x) —a, sin x — a, sin 3x . . . —a,sinnx)? dx . (15) 
z Jo 
Differentiating, we see that this is a minimum by the variation of a, when 
a/2 
= (f (x) —a, sinx —a, sin3x .. . —a,sinnx)sinnxdx (16) 


This expression represents the ordinary rule for determining the coefficients in 


| 
| 
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Fourier’s series by multiplying by sin na and integrating. Moreover, the mean 


square of the error is equal to 


== { f (x) } — 4 (a,? (1 


a 


= 2 4 q2 


We have chosen a series of odd sines as an example, but the same method applies. 


to any series. From this the degree of accuracy of any approximation can be 
estimated. 


It is important to notice that according to the method of least squares, the 
sum of any number of terms of the series, or even a single term or a collection 
of terms picked out arbitrarily gives the best possible approximation to the 
function when the values of the coefficients are those determined by Fourier’s 
Theorem. For example, the value of P, sina will make the mean square of the 
error a minimum in the interval from o to 90° when it is sought to replace f (a) 
by a single term containing sina. It therefore gives the best approximation for 
the whole range of values from o to go°. It is, therefore, not however the best 
value to take when only small angles have to be considered. 


But the method of least squares can also be used to obtain the series which 
gives the best approximation to a function f (a) over any interval whatever. The 
coefficients will still be obtainable by the ordinary method of multiplying and 
integrating as is evident from equation A. The only difference is that unless the 
limits of integration are o and 90°, or similar angles, the integrated products of 
sines or cosines will not vanish and the coefficients will therefore have to be found 
by solving a number of simultaneous equations. Examples of this method will 
be given later on. The following application may be found useful in aero- 
dynamical calculations. 


PrRoBLEM.—Given the value of f(a) in the form of a Fourier’s series between 
the limits o and 90°, to find the value of A in order that the single term A sina 
should give the best possible approximation to f (a) in a given limited interval, 
say, between o and 30°. 


Taking the original given expansion to be XP, sinna with n odd, and 
assuming A = P, + K, the method of least squares gives 


sin?’x — P, sinx sin3x— P,sinxsin5x ...)dx =o. (18) 


For the limits o and 30°, which are often required in aeronautical work, | 
find that 
3 \ P, P, P, P, P,; 
K | ———1] = — + — + ——— — — — (19) 
9 ] 2 2 4 20 5 


It will probably be better to employ this or some similar formula where it 
is required to employ the sine rule for small angles of attack rather than to use a 
series, the convergence of which is doubtful, for the limit of f(a)/sina when a 
approaches zero. 

16. Formule for the Centre of Pressure.—According to ‘‘ Joessel’s Formula ”’ 
the distance of the centre of pressure of a plane lamina in front of the centre of 
the area is given by 

0.3 (I — sin a) 
referred to a lamina of which the whole breadth is unity. In some books the 
coefficient is quoted as 0.33 instead of 0.3. 


— 
| 
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This formula can only hold good if a is in the first quadrant. If a 
is in the second quadrant it is obvious that the corresponding form must be 


—o.3 (I — sina). 


To obtain a formula which shall hold good in all four quadrants it is easily 
seen from considerations of symmetry that the expression must be expanded in 
cosines of odd multiples of a between the limits o and tz. We obtain on multi- 


plying and integrating 


} (cosa COS 3a COS 54 COS 7a 
0.3 (I — sina) = 0.3 » | —— 
\ 233 5-6 6.7 
COS Oa COS Ila \ 
9.10 10.11 ] 


According to the Kirchoff-Rayleigh theory of discontinuous liquid motion the 
corresponding expression is 


and the coefficients in the Fourier expansion may be determined by a formula of 
reduction similar to that described in £ 8. 
Another formula, sometimes described as Soreau’s formula, is 


I 


4(1 + 2 tana) 
which may be written 


4 (cosa + 2 sina) 
Che integrals required for evaluating the Fourier coefficients may be evaluated 
by the reduction formula 


5 cos na — 6 Cos (n — 2) a + 5 cos (n — 4) a 
———— = 4 cos (n — 2) a (cos a — 2 Sin a) 


cosa + 2Ssina 


(21) 


Che numerical values of the coefficients in these expressions have been formed 
by Mr. Creak and are given in Table G. 


[he corresponding coefficients calculated from several of Eiffel’s experimental 
results for plane areas are given in the same table and have been calculated by 
Mr. Caradog Williams. The method emploved in this case is the same as in 
the calculation of the pressure coefficients. If @ (a) is the function to be expanded 
and we assume 


y = (a) = C, cosa + C, cos 3a + C, cos Sac . 
we obtain 
(x) sinnx (x) sin nx 
+ n jo n 
Jo fe) 


‘ ‘ (22) 


In this case @ (7/2) = o so that the terms in square brackets therefore vanish. 


| 
COS a 3 
4 
I 4 
cos a 
| 
4 
L 
4 
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If now in the interval between two observed data we replace @! (x) by its 
mean value Ay/Ax, we obtain 
[ Ay 


cos NNy+1 — COS NX, 
>, 


4 | Ax 


J I Ay \ 
SNcosnx A | -—- 
\ Ax | 


In this formula, angles are measured in radians. In practical work Aw is 
necessarily measured in degrees, so that when thus measured the formula 
becomes 

4 x 180 Ay |r 
C - nx,4 | —— — |—— (23) 


A difficulty arises from the fact that when a o the pressure vanishes and 


| 
(| 
consequently @ (0) is indeterminate. Consequently the lowest value of a for which 
@ (a) is tabulated in Eiffel’s experiments is 2°. This difficulty is, however, not 
so serious a one as might have been expected. The only indeterminate quantity 
is the term 
COs 2° - [ 


Unless n is large, cos n 2° I is very small, and so this term is not very 
important. Mr. Caradog Williams has dealt with this question in two different 
ways, first, by neglecting this term altogether, and secondly, by assuming: that 
the value of Ay/Ax between the first two tabulated data remains the same down 
to a o. In the former case the curve is assumed to be horizontal from a O 
up to the first tabulated value; in the latter it is replaced by straight line from 
a o up to the second tabulated value. 


The difference in the numerical values of the coz@ffcients is however small. 
For example, in the case of the rectangle 15 x 45cm., aspect ratio 4, we have 
> 


33 and that in the second case 


found that in the former case C .009 
C,,= .007783. The expansions converge rather more slowly in the case of the 
experimental data than in the theoretical formula. This is, of course, due to the 
fact that according to Eiffel’s curves @!' (a) is considerable when a is small and 
of course it changes sign with a. The formula may therefore be expected to 
agree fairly well with the assumed tabulated data except when a is small. 
17. Formula for the Moment of Resultant Pressure.—Assuming Duchemin’s 
formula for the resultant pressure and Joessel’s formula for the centre of pressure, | 
the moment of the resultant pressure about the centre of figure is given by 
2 2 sina 
- (1 sin a) 5 (24) 


10 I sin? a 


| 


4», Would have when 
acting on an arm whole breadth of plate. | 


taking as a unit the moment which a perpendicular thrust P 


As in $4 the moment of the resultant pressure can be expanded in an even 
ime series | 
Sin 2na 


where 


\ 


| | 
- | (1 — sina) sin 2 nada 


10 + be 
re) 


XUM 


} 
| Ax | n? 
| 
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The formula of reduction in this case becomes 
sin na —-6 sin (n — 2) a + sin (n — 4) a 
=—4sin(n—2)a. (25) 
I + sin? a 


The numerical coefficients in the expansion have been calculated by Mr. Creak 
and are given in Table H. 

As in previous examples, however, if approximate values are taken the 
coefficients calculated from this reduction formula ultimately form a divergent 
sequence. 

The calculation of the corresponding series for other theoretical formule, as 
also for the results of experimental data, presents no special difficulty and it 
hardly appears worth while to consider any further cases in detail. 

18. In a subsequent paper we propose to discuss the best way of applying 
Fourier’s series to the aerodynamical coefficients of cambered planes, aeroplane 
wings and complete systems of aeroplane surfaces. Such cases present a number 
of peculiarities which will require careful consideration. 

It would be premature to express any opinion as to the practical value of 
the Fourier method, because everything must necessarily depend on the class 
of investigations for which their utility may have to be considered. All that the 
present paper claims to do is to give a general idea of the kind of results to which 
this method leads, and any special peculiarities presented by the methods most 
convenient for obtaining them. It will, I think, be seen that the methods here 
described involve several points of difference from the subject matter contained in 
existing papers and text books. It may be further stated that most of these 
differences were not anticipated when the present investigation was commenced, 
with the result that what was thought to be a short piece of work has grown 
into a long one. 


TABLE A.—PRESSURE COEFFICIENTS (THEORETICAL). 


Pali, - SP, sin na P P. P. P. P | P 

Duchemin’s formula 1°1716 ‘2010 "0345 "0059 ‘OO10 ‘0002 
Gerlach’s formula | 

(modification of 

Rayleigh) 1'0794 ‘0979 0252 ‘0097 ‘0047 "0926 
Renard’s formula 

(modified for P,,=1), 1°2893 "1357 \—'0742 ‘0208 —‘o129 ‘0078 
Soreau’s formula "25 
Eiffel’s formula (1907) 

Poa up to 30° and 

P=1 from a= 30° to 

a=90° "2702 ‘0486 —'0248 —'0300 
P=2sina up to 30° 
P=1 from a= 30° to 

a= 90° ‘2757 0551 —'0197 —'0276 |—‘oI00 


he 
| 
> 4 

| 
i 
| 
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TABLE B.—PRESSURE COEFFICIENTS. 
EIFFEL’S EXPERIMENTAL RESULTS FOR PLANE AREAS. 


Pa/P,,=>P, sin na 


(n odd). 


P P, p P. 


Square Plate 25 x 25cm. 


Aspect Ratio I 173103 ‘4137 ‘0794 —'OdSI5 —'0830 —0024 
Rectangle 22.5 x 15cm. 

Aspect Ratio £0554 "2055 "1385 0372 —'0208 —'04y3 
Rectangle 30 x 15cm. 

Aspect Ratio 2 ..| ‘1299 "1342 "0313 -—'009g0 
Rectangle 

> 

Aspect Ratio 3 . 10059 "2193 "1157 "0572 0335 0058 
Rectangle 90 x 15cm. 

Aspect Ratio 6 | *2000 1264 ‘O741 0478 0254 
Rectang k go x 1ocm. 

Aspect Ratio 1°1462 "2399 "131g ‘0518 ‘0309 
Rectangle 15 x 45cm. 

Aspect Ratio 4 1°2001 "25600 —'0733 —'0893 —‘oo8o0 ‘O214 
Rectangle 15 x gocm. e 
Aspect Ratio é 1°1749 ‘1439 —'0878 |—-0361 

TABLE C. 
n I 2 5 7 Q 11 13 15 17 IQ 21 
} x 180 

leg ___| 1°8631 “9089 1652 1729 | 1°9546 1°7803 1°6352 | 1°5109 | 1°4022 1°3056 1°2187 
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TABLE D.—SPECIMEN CALCULATION. 


SQUARE PLATE 25 x 2icm., ASPECT RATIO 1. 
a Yn Quotient. Factor F. log F. 
oO | —= 
"105 5 033 
5 o'105 | 3°7702 
5 °039 
10 0°36 | | —"0055 3°7404 
10 
20 o'805 | OO10 3/0000 
“435 10 0435 | 
: = 
30 1°24 OLI5 2°0007 
| 10 5 
35 1°40 | ‘O12 2°0792 
(ole) 3 O2 
» . | . ~~ 
35 1°46 105 2075 
2 —145 
40 —'127 1103! 
5 —o18 
45 1°08 —'O147 2°1673 
—'05 15 —'0033 
60 1°03 —0023 3°3017 
15 —'oo! 
75 1015 ‘000 
O15 rs OO ais 
go | Jj 3°0000 
CALCULATION OF P.,. 
a 5 10 20 30 35 | 38 40 45 60 go 
3a 15 30 60 90 105 114 120 135 180 270 
Nearest acute 75 66 60 5 oO j.90 
omit 
log sin 30 14130 1°6990 1°9375| 1'Q849 1*Q007 1°9375 (logo) 070000 
log F 3°7782| 3°7404| 3°0000| 2°0607] 2°0792] 1°2175| j 1°1038| j 2°1673 | j j 3°0000 
log X j 2°1001 | j 2°3483} 3°8464] 2°9696| 2°9730) j 1°9502) j 1°9257 9 
0'01259 | -0°0223 |0°007021 | 0°09324 | 0°09397 1°222, -0°8915 |-0.08428 0°008 108 


Positive ‘007021 Negative — ‘01259 
“093240 2230 
‘093970 — ‘8Q150 


*222000 


008108 


XUM 


424330 


‘010070 


"413669 = 


1007 


g 
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TABLE E.—PRESSURE COEFFICIENTS (COSINE SERIES). 


Pa/P,,=sin aXQ.,, cos anal Q, Q, Q, | Q, 
Square Plate 25 x 25cm. | 
Aspect Ratio I | 1°6703 °6883 | —"1354 | —'2245 | —'0833 0404 
Rectangle 22.5 x 15cm 
Aspect Ratio | "2016 | —*0682 —-1450 | "0986 
Values of P,,,, obtained from above. (Compare Table B.) 
| 
; 
I n+1° coefiicient of P P. P. P. P. P,, 
sin (n + 1)a. 
Square Plate 25 x 25cm. 
Aspect Ratio I 1°3301 “ATIQ “O445 ‘0700 —0649 — ‘0086 
Rectangle 22.5 x 15cm. 
Aspect Ratio = 1.5 [1774 2015 "1349 | —'0232 —"0440 
TABLE F.—LIFT COEFFICIENTS. 

Square Plate 25 x 25cm. 
Aspect Ratio I ‘8650 2405 —OOII |—'o0822 —'0427 
Rectangle x t1ocm.! 
Aspect Ratio 9 6030 "1859 "1054 ‘0054 0413 279 
Rectangle 15 x 45cm. 
Aspect Ratio 1 ... 7604 |» ‘0907 —0813 | —'0487 ‘0007 ‘O10! 

DRAG COEFFICIENTS. 

Drag /P,, = cos 2na. D, D. dD, 
Square Plate 25 x 25cm. | 
Aspect Ratio I “6582 |—"4513 —"1671 | — 0805 | —-0007 "0403 
Rectangle go x 10cm. 
Aspect Ratio 9 "4531 | — 0540 |—'0264 | —'0136 | —‘o104 
Rectangle 15 x 45cm. 
Aspect Ratio ...| °6330 ‘5057 | —0907 ‘or80 *0407 ‘O147 


XUM 
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TABLE G.—CENTRE OF PRESSURE COEFFICIENTS. 


(Distance of centre of pressure from centre of figure 


full breadth of 


plate x XC, cos na (n odd). 
| 
C, C, C; C, C, 
Joessel’s formula “IQIO "0637 “O27 ‘0042 0035 
Rayleigh’s formula "1439 ‘0200 ‘0075 0038 0023 .OOI4 
Soreau’s formula 1555 O3I1 0176 ‘0094 ‘0000 ‘OO44 


EIFFEL’S EXPERIMENTAL RESULTS FOR PLANE AREAS. 


Square Plate 25 x 25cm. 


Aspect Ratio i 1790 
Rectangle 45 x 15cm. 
Aspect Ratio 2 | 
Rectangle = 90 x 15cm. 
Aspect Ratio 6 388 
Rectangle 15 x 45cin. 
Aspect Ratio 1 ..| “2054 
Rectangle 15 x gocm. 
Aspect Ratio 1 | "O26 


0553 ‘O401 
0301 ‘0284 
‘0261 "0235 
‘OO1O 
—'0043, —'0055 


‘O145 


‘OIIO 


‘0097 "0052 
O152 "0093 
‘0126 ‘0089 
0240 ‘0078 
0082 —'0035 


TABLE H.—MOMENT OF RESULTANT PRESSURE COEFFICIENTS. 


(Moment of resultant pressure about centre of figure 


unit the moment of the normal pr 
bread 


ressure P,, 


th of the lamina.) 


XG,,, sin 2na taking as 


acting on an arm equal to the full 


G, 


From Duchemin’s and 
Joessel’s formulze ...) *0935 
From Soreau’s formule | 
| 


"0950 


Ga Gi, 
0515 ‘OI1g | *0055 
0267 | ‘0082 003 


‘OO16 


‘O002T 


G,. G,, 
‘0013 ‘0006 
‘OOII 


‘0225 
' 
‘0173 
‘O131 
= 
XUM 
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LONG TAILS VERSUS SHORT TAILS IN 
AERODYNAMICS. 
RY PROF. G. H. BRYAN, SC.D., F.R.S., (BANGOR), 
AND 


S. BRODETSKY, M.A., PH.D., (BRISTOL). 


1. Consider any form of aircraft provided with a small rudder plane placed 
behind its centre of gravity and at a certain distance | from it, which we may 
speak of as the length of the tail. When this tail is subjected to a given wind 
pressure, the moment of this pressure about the centre of gravity, tending to 
turn the machine round, is proportional to 1. From this point of view a long 
tail is preferable to a short one. 


On the other hand, the pressure of the air on the tail plane depends on the 
normal velocity of the wind relative to it. Hence the normal velocity of the plane 
itself is limited in magnitude, and the angular velocity with which the machine 
turns about its centre of gravity must therefore be less when / is large than when 
l is small; that is, the long-tailed machine turns more slowly than the one with 
a short tail, and this is a disadvantage. 


We thus have two conditions which are to a certain extent incompatible, 
and it becomes interesting and useful to work out a problem in which the joint 
effect of these conditions is taken into account. In order to avoid complications it 
appeared to be desirable in the first instance to consider the simplest possible form 
of system furnished with a tail. The system discussed in the following paper might 
perhaps be described as a ** dart,’? and may be roughly illustrated by an ordinary 
arrow. It is assumed to be a heavy body experiencing no resistance itself from 
the air, and acted on by no forces, but furnished with a rudder plane, the 
air pressure on which follows the sine law of resistance when the plane is moving 
in a direction other than tangential relative to the air. In this case the pressure 
on the tail is wholly employed in overcoming the inertia of the system ana 
changing its motion so as ultimately to bring the direction of motion and the 
line joining the plane to the centre of gravity into coincidence; or, in ordinary 
language, to make the system fly in the direction in which it is pointing. The 
problem presented by an aeroplane or an airship is of course much more com- 
plicated, because the pressure on the rudder is largely employed in overcoming 
the air resistances on the machine itself. In the case of an aeroplane which cannot 
turn curves without banking, the discussion of the action of a vertical rudder 
brings us back to the problem of lateral oscillations and lateral stability, while 
the effect of a horizontal rudder is considered in connection with longitudinal 
stability. 


2. A plane lamina, centre C, so narrow that for moderate angles of attack 
we may use the sine law for the air pressure and ignore the shift of the centre of 
pressure, is attached to a system so that the whole mass is .\/, and the centre of 
mass is at G, outside the lamina but in its plane. Originally the system is moving 
with a given velocity V along the lamina and with no rotation. A wind of 
constant velocity IV is suddenly generated at a small angle a with the direction of 
motion, backwards. To find the oscillation set up, and the best position for the 
tail, whose size is given. ‘The motion is supposed two dimensional—in the plane 


of the paper—-and there are no applied forces. 


It will simplify the work if we consider, not the actual motion of the system, 
but the motion relative to the air. The initial relative velocity can be found by 


XUM 
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the following construction (see Fig.). Measure off CK containing W units of 
length, so that KC represents I, the velocity of the wind; and draw KB parallel 
to the initial direction of CG containing V units of length, and therefore repre- 
senting the initial velocity V of the dart. Then CB represents the initial velocity 
U of the dart relative to the air, and similarly BC, which is equal and opposite 
to UB,* represents the relative velocity with which the air appears to blow on 


G B 


R 


the lamina initially. The problem of the relative motion thus becomes the same 
as that of a dart originally projected in the direction (B with a velocity U 
measured by the length CB. If wv are the relative velocity components of the 
dart along and perpendicular to the lamina, 6 the inclination of the lamina to the 
fixed direction (b, 6, the initial value of 6 so that 6, angle ABC, then we 


have initially 


uw, = Ucos6, = Wecosa + JV, 
v, = Usin 6, = W sina; 
and if a is small, we may take 
W 
V; 6. — 25 05 = Wa . 


Che quantities u, v, 6, define the motion of the dart relative to the air. It 
is clear, however, that 6 represents the actual direction of CG at any moment. 


3. If CG, the distance of the centre of the lamina behind the centre of mass 
of the system is 1, and Mk? is the moment of inertia of the whole system about 
an axis through G normal to the plane of the motion, the equations of motion 
ere as follows :— 


du dé 
dv dé d?6 R 
dt dt dt? M' 


where Ff is the air resistance. As the lamina is narrow, we assume the sine law 
of resistance and put 

R 

M 


Av (uw? + v?) 1/2, (3) 


Note.—If the velocities of two moving bodies be represented in magnitude and direc- 
tion by two sides of a triangle, both drawn from their common vertex, then the third side repre- 
sents their relative velovity. Thus in the above figure where KC and KB represent the velocities 
of the wind and dart, the relative velocity is represented by CB or BC according to whether 
we are considering the relative motion of the dart or of the wind. This simple construction for 
relative velocity can be easily verified by supposing two bodies to start simultaneously from 
the point K with the velocities Wo and V; then in a unit of. time they will have arrived at 
C and B, so that the line CB will show the change of relative position in unit time, and will 
therefore represent the relative velocity. This explanation appears necessary because the treat- 
ment of relative velocity in many text books is very confusing, and it is our experience that 
even advanced students frequently make mistakes in consequence. 


Vv 
K 
V 
‘6 
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where A is the reciprocal of a length depending upon the size, shape and aspect 
of the lamina. 

We suppose a small, so that 6 remains small for a finite oscillation, and v 
is also small. We therefore neglect squares and products of d6/dt and v, and 
the equations (2) become :— 


du 
dt 
Al 
3 Uv 
dt? ke 


This means that w is constant to the first order of small quantities and equal 
to U. Eliminating v we get 


d°6 AU (I? + k?) d?6 AlU? dé 
+ — 4 — (5) 

dt* k? dt? k? dt 

The solution of (5) is 
— 
where ,, A,, are the solutions of the algebraic equation 
AU (I? + k?) AlU? 

k? k? ( 


and P, Q,, Q., are constants defined by the initial conditions. 


Measuring the time from the beginning of the disturbance caused by the 


wind W, we have when t = o, 
dé 
6 = =o, v= Usin#, U6,, 
dt 
so that at t = o, 
d*6 AlU? 
— = — = : (8) 
dt? k? 
to our degree of approximation. After some reduction we find 
A,t — 
A.€ 6 (o) 
) 


It follows from (7) that for a stable oscillation we must have | positive, and 
when this condition is satisfied, i.e., when the lamina is behind the centre of 
mass of the system, the motion is always a stable one, the system gradually 
attaining some different steady motion. The ultimate direction of relative motion 
is along CB. The effect of the tail in this case is to alter the orientation of the 
dart without affecting its relative velocity. The ultimate orientation is also along 
CB, t.e., along the direction of the initial relative velocity. 


4. The solutions of (7) are 
AU (I? + k?) A?U? (l? + k?)? 
As + | 4k ke | (10) 
The advantage of a very small value of !| is that no actual oscillation arises, 
but a gradual change of direction. But in such a case the decrement of the 


XUM 
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deviation from steady motion would be slow. In practice A,, A,, would be complex, 


so that a damped oscillation ensues. If we write 


Nis A. = Pp iq; 
where 
AU (I? + k?) 
— 
and 
AlU? 12U? (1? + k?)?) 
q 4k* j 
we obtain for (9) 
6= 6, ; sin ( qt + tan ) (12) 
Thus at any moment the amplitude of the oscillation is proportional to 
4 2\1/2 


q 

In order to get the quickest decrease in the amplitude, no attention being 
paid to the initial extent of the oscillation, we must make p as large as possible, 
i.e., by (11), we must make /?//° as large as possible. If, on the other hand, 
the important thing is to avoid violent oscillations, then we must make 
(p? + q*)'/?/q as small as possible. This condition easily reduces to (I? + k?)?/k?l 
as small as possible. 
5. In an ideal case, where the tail and the connecting rods are supposed to 
have no mass, k? is independent of |. We then get that in order that the 
oscillation may die down soon, | must be as large as possible. But in order to 
avoid violent oscillations, we easily find that | must be as nearly as possible equal to 
I 3. For practical purposes we conclude that | must be at least equal to k/./3 


A 
since diminution below this limit involves a double disadvantage. 


In an actual flying system, we may expect that 
k? =a + : (14) 
where a and b are known positive quantities. We again get that to make the 
oscillations die down soon, we must make | as large as possible. In order to 
prevent violent oscillation, we must have 
(l? + a + DI?)? 


TS 
(a + 1 (15) 
a minimum. This expression is stationary for values of | given by 
[2 
b(r + +3— —1=0, 
giving the unique possible value 
—3+(9 + 4b + 4b?) 
a 2b (1 
Thus the oscillation is least violent for 
} 
2b (1 + b) J 
since this makes (15) a minimum. In practice this value of | would be a lower 
limit. ° 
When b/a is small, (16) becomes approximately 
Ja 


a/ 
agreeing with the result previously found, 
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6. Since 6 ultimately vanishes it will be seen that the direction of motion ulti- 
mately becomes parallel to CB, and this therefore represents the direction of the 
asymptote of the path. It is, however, important to determine the distance of 
this asymptote from (, as this represents the total displacement of the dart in 
a direction perpendicular to ('B relative to the air. This displacement is obviously 
equal to 


sin 6— v cos 6) dt é (18) 


which for our purposes of approximation may be taken to be 


(U6 —v) dt 
Now from equations (4) we obtain 
dv dé + 12 
dt lt (19) 


This is the equation obtainable by taking moments about the centre of the 
tail. It follows from this that the total displacement is equal to 
— 0, + Ct (20) 
where (, and C, are constants. But since the initial velocity perpendicular to 
CB is zero and 6 is initially equal to @,, the total displacement becomes equal to 


3: 72 (21) 
l 
This is a minimum when [* = k?. If, however, k? a + bl*, the condition 
for a minimum becomes a = bl? + I? and therefore k? (1 + 2b)I?. In any 


case there is a certain length of tail for which the displacement is least. 


If, however, we consider the total displacement of the centre of gravity we 
must subtract 16 from the above amount, and the total displacement sideways 
now becomes equal to 6,k?/1, showing that a long tail is the best. 


It is easy to obtain a generalisation of this last result for the case of any 
kind of aircraft whose motion may be treated as two dimensional. In such a 
case there is at every instant a point on the axis of the system through which 
the resultant air resistance passes. This point takes the place of ( the centre 
of the tail. The main difference is that since the components of air resistance 
depend both on the linear and angular velocities of the system, and these are 
not in general proportional throughout the motion (for v is initially finite and 
d@/dt initially zero), the distance of this point from G, which we have hitherto 


denoted by | will now be variable. The total displacement of G at right angles 
to the direction of final motion will now be 


We may write this equal to k76, {l] where [1] denotes some quantity such that 
1/[{l] is intermediate between the greatest and least values of 1/I. In any case 
it is advantageous from this point of view that [/] should be large. 


7- The preceding investigation of the small oscillations is, however, only a 
particular case of the more general problem of aeroplane stability. A more useful 
estimgte of the initial action of the tail plane can be obtained by expanding @ in 
powers of the time by means of the exponential series. We thus obtain from (9) 


>! > 1 
2. 3 


4! 


7 
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Making use of the formule for the sum and product of the roots of the 
AU (1? + k?) 
4 

(24) 


24 


quadratic (7) this gives 
ALU? 

6,—6 é 

hk? 

AW LAE 

This equation shows that the initial angular acceleration of the system is 
greatest when /| is large, but in this case the angular acceleration decreases more 
rapidly than it would do if 1 were small, thus confirming the general statements 


at the commencement of the present paper. 
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ON THE FLIGHT OF LOCUSTS. 
BY DR. E. H. HANKIN, M.A., SC.D., A.F.AE.S., AGRA, INDIA. 


Some time ago when leaving my laboratory I happened to notice a large 


number of cheels gliding overhead that seemed to me to show something unusual 
in their mode of flight. I guessed that they were after locusts, a flight of these 


insects having passed over Agra three days previously. 


On reaching my house 


I found, with the help of my binocular, that my conjecture was correct. The 
following are the notes that I took at the time: 
1ith August, 1915, 3.0 p.m. A few small thin clouds. Air rather hazy 


at a height. Wind S, occasionally moving small braaches. 

A large number of cheels were gliding in a way that seemed to me 
unusual. They were at heights between 200 and 400 metres. They were 
generally slow-flex-gliding and did so in all directions and with frequently 


changing course. Two were seen rapidly losing height, feet foremost. 
Occasionally they dived or at least glided steeply downwards at speed. I 


guessed they were hunting locusts and with my binocular soon saw _ that 
several of the cheels kept bending down their heads to eat locusts that they 
held in their claws. 

3-6. A locust noticed in flight. It was about 100 metres overhead. 
Flapping was with all four wings. It was flapping to the west, but was 
being carried almost due north by the wind. 


z 


Three locusts seen in flapping flight. 


3°7° 
3-10. <A locust seen flapping at a height of about 100 metres. It stopped 
flapping and began gliding with loss of height. Its wings were dihedrally 


up (The loss of height was very well marked. It may have been at the 
rate of 1 in 5.) 

3-11. A cheel seen gliding steeply downwards for a distance of about 
50 metres. Then it caught a locust in its claws and at once changed te 
slow speed horizontal flight. Making a measurement with my telemeter I 
found it to be at a height of 300 metres. I could see the vellow colour ot 
the locust while held in the claws of the cheel. 

Besides the above noted case I saw two other locusts gliding with loss 
of height. In one of these cases the wings were clearly seen to be dihedrally 
up. Otherwise the locusts were, as usual, in continued flapping flight. On 
the other hand, the cheels, whether soaring with gain or loss of speed or 
with gain or loss of height, were always gliding in their efforts to catch the 
locusts. They aimed directly at the insect and caught it as usual in then 
claws. Occasionally a slight irregularity of flight was noticed just before 
the cheel caught the locust, but in such cases the cheel was seen to make a 
sudden change of course as if it had only just seen its prey. 


The locusts were not in a thick flight. I never noticed more than two 
in sight at once. The cluster of cheels kept moving their position, evidently 


being able to see the locusts from several hundred metres distance. 


A few days later I made a similar observation. It was as follows: 


16th September, 1915, 8.45 a.m. In Hewett Park, Agra. A few 
locusts about. Cheels were catching them at heights between 4o and 150 
metres. They did so while in continued gliding flight. On two occasions I 
saw the cheel make a sudden rotation up round its transverse axis at the 
moment of catching the locust and occasionally they showed a single flap, 
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and once three flaps, after catching the insect, as if to regain balance or 
possibly speed. In gliding down from a height to the insect the cheels always 
travelled in an apparently straight line. One showed a slight lateral wobble 
at the moment of seizing its prey. Thus, their flight, with the above excep- 
tions, was entirely gliding whether they were gaining or losing height. 


Soaring and Non-Soaring Flight at the Same Time and 


Place. 

Obviously the interesting point about the above observations is that the 
soaring bird and the non-soaring locust were flying together in the same air. 
The feebly flying locusts were all being drifted to the north despite their industrious 
apping in a westerly direction. In every case their course appeared to be a 
straight line. At the height above the earth at which they were flying in a light 
wind there were no local currents capable of affecting their course to an appre- 
ciable degree. Nevertheless, the apparently homogeneous air in which they flew 
was well suited to the soaring flight of the cheels. In the same air one flying 
creature could only maintain itself by constant flapping, and at once lost height 
if it glided, while the other maintained its speed, which was a greater speed, by 
simple gliding whether this gliding was accompanied with gain or loss of height. 
In a similarly light wind at tree-top level there are probably ascending currents 
and other air irregularities capable of attecting flight. At this level cheels, when 
catching flying white ants, always use a great deal of flapping in their flight. 


As it does not often happen that soaring and non-soaring fliers are seen in 
flight side by side the following instances are worth quoting :— 

2oth August, 1915, at 11.30. Sunshine. 

A group of cheels and one eagle were seen circling at a height of about 
150 metres. Two crows were present and were teasing the eagle. They 
were in continued flapping flight. The eagle responded chiefly by a sort ol 
slow double dip movement of its wings. I watched this for several minutes. 
The wind was west and very light as while I was observing the group of 
birds did not drift very far away from overhead. At length the crows dived, 
with wings flexed and retracted, coming down in a swerving course. The 
crows had been flapping till the moment they dived while the eagle and the 
cheels were in continued gliding flight. While this was going on two soaring 
birds of unknown species joined the group. Their wings appeared snow 
white except that the anterior margins appeared grevish-vellow. A vulture 
passed flex-gliding just over them at about 200 metres height with the white 
part of its wings appearing bright vellow. 


In this case the crows were above the earth at a height that was very unusual 
for them. This was not the case in the following observation :— 

19th July, 1915. Bombay. 8.0 a.m. Slight sunshine. Three crows 
seen in flapping flight and one cheel slow-flex-gliding. They were all going 
in the same direction and at a height of about 10 metres above the house tops. 
The wind was light and on their beam. Thus, the same air (the cheel was 
closely following the crows) supported one bird in soaring flight while the 
other birds had to flap. 


In this latter case the non-soaring birds travelled at the same rate as the 
soaring bird. In the former case they appeared to travel faster. The converse 
relation may be observed. For instance, I have seen a crow teasing a cheel at 
the time of commencement of morning soarability. The cheel was circling with 
an occasional flap. After a few minutes as the air became more soarable the 
cheel circled without flapping and then glided away at a speed with which the 
flapping crow was unable to keep up. Some time after the air has become soar- 
able for cheels it becomes soarable for birds of heavier loading such as vultures. 


= 
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While this is happening vultures are liable to be teased by cheels. The following 
is an instance :— 
11th October, 1913. Calcutta. At 9.55. Air thick. Isolated cumulus. 

Sunshine. Wind north, slightly moving leaves. A black vulture at 300 
metres height being teased by a cheel. The vulture appeared to have a 
slightly greater speed. It was circling (without flapping) most of the time. 
The cheel kept approaching by cutting off corners and describing smaller 
circles. When near the vulture the cheel usually made one or two flaps. 
The black vulture responded by dipping movements of both wing-tips. 
Occasionally it made a large momentary canting resembling a strong depres- 
sion of one wing. On reaching a height of about 500 metres the black vulture 
slow-flex-glided away, being chased for about half a mile by the cheel, which 
it slowly outdistanced. The teasing had gone on for about 10 minutes. 
Thus, in this case, when the air was fully soarable the bird of heavy loading 

was able to glide horizontally at a greater speed than the bird of light loading. 


Soaring and Non-ssoaring Dragon-flies. 


But perhaps the most interesting case of soaring and non-soaring animals 
flying in the same air is afforded by the habits of the soaring dragon fly, Pantala 
flavescens, which in Calcutta often flies in company with another species of 


dragon fly, Rhyothemis variegata. This latter species has a restricted power of 
soaring, but is usually seen in flapping flight. These flapping dragon flies fre- 
quently fly together in a group consisting of from a few dozen to a few hundred 
individuals. The group remains for hours together in the same spot, often to 
leeward of a large tree. The individuals in the group keep rather close together 


so that any one dragon fly is rarely more than two or three feet distant from its 
nearest neighbours. Occasionally a few specimens of the soaring dragon fly 
may be seen gliding in the group. The following description is quoted from my 
book, ‘‘ Animal Flight,’’ page 388 :— 

“It is noticeable that the Pantalas when gliding in a group with 
Rhyothemis travel at a much greater speed than the latter. That is to say, 
in the same air and within a few inches of one another, one dragon fly uses 
flapping in its flight. The other dispenses with flapping, but nevertheless 
travels at perhaps double the speed. An upholder of the ascending current 
theory .might devote -much learned—and perhaps useless—casuistry to ex- 
plaining why the swift-gliding Pantala finds these supposed ascending currents, 
while the flapping Rhyothemis so consistently avoids them. The facts of the 
case contradict any theory that seeks to explain the gliding of Pantala by 
means of pre-existing movements in the air. Supposing that Pantala gets 
energy from the air by meeting eddies, or from some mysterious aerial vibra- 
tion, or from some turbulent movement of the air, why is it that Rhyothemis 
is unaffected thereby though a lighter insect? ”’ 

In view of the above observations it is obviously desirable to accumulate data 
that will enable us to compare the flying capacities of different kinds of flying 
animals. 


The Flight-Speed of the Locust. 


With regard to the speed of the locust, observations are very easily made. 
A locust when liberated usually flies away horizontally for a considerable distance 
in a straight line. It is therefore easy to time it over a measured course. I found 
in this way that the speed of flight of locusts, both in my laboratory and in still 
air out of doors in the early morning, is 4 metres per second. 


Accurate data are yet lacking for making a comparison between the speed 
of the locust and that of other insects. As compared with most butterflies its 
flight appears to be slow. In one experiment 1 found that a soaring dragon fly— 
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Tramea burmeisteri—could fly upwards by flapping on a course making an angle 
of about 60 degrees with the horizon at 4 metres per second and that when it 
had reached a height of about 3 metres it glided horizontally without flapping at 
a distinctly greater speed. 


The Wing Loading of the Locust and Other Insects. 


The following table permits of a comparison between the loading of the 
locust and that of other insects :— 
TABLE I. LOADING OF INSECTS. 


LOADING, 


ORDER, NAME. LBS. PER KILO. PER 


SQ. FOOT. SQ. METRE. AUTHORITY. 


ORTHOPTERA. 


Locust .| 125 | 61 | H 
Dragon flies— | 
Rhyothemis variegata ‘| | H 
Maiden dragon fly [ 029 | 
Pantala flavescens | 16 | H 
Small dragon fly 16 | del 
Common dragon fly 046 | 23 | 
Hemianax ephippiger cen 066 | .32 | H 
Flat-bellied dragon fly sa 094 | .46 | LE 
| | | 
NEUROPTERA. | | | 
Ant lion H 
| | | 
DIPTERA. | | | 
Gnat | 020 | | L 
Daddy long legs (Tipula) O08: | | 
House fly | 14 | .56 | 
| | | 
LEPIDOPTERA. | | | 
Cabbage butterfly | wer | - | 
Argynnis childreni | | «ta H 
Papilio ravana | 028 =| | H 
Swallow tail butterfly... 063 | L 
Hawk moth (Cherocampa celerio) | 67 | & | H 
Sphinx moth | 214 | 1.05 | E 
| | | 
COLEOPTERA. | | | 
Cockchafer | 195 | 96 | I 
Stag beetle (Lucanus cervus), female | 21a, | Ios I 
male | 266 1.30 | 
Rhinoceros beetle 318 1.56 I 
Sternocera_nitidicollis, calculated on all four | | | 
wing's | 3.60 H 
hind wings only 1.34 6.57 H 
| | 
HYMENOPTERA. | | | 


Note.—Measurements marked L are from Lougheed (‘f Vehicles of the <Air,’’ p. 162). 
I 

Those marked de L are based on data given by de Lucy quoted in the article on Flight in 

the ‘‘ Encyclopedia Britannica.’"”. My own measurements are marked H. 
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An interesting point to notice in this list is that the beetle Sternocera has a 
loading comparable to that of a vulture. When in flight this beetle makes a 
remarkably loud humming noise. Its speed ahead is very slow. Hence, it is 
improbable that the fore wings or wing cases aid appreciably in support. Another 
reason for holding this view is that in flight they are held at a dihedrally up 
angle of about 45°. Consequently it is probably correct to regard the figure 1.34, 
which is calculated on the area of the hind wings alone as the actual loading. 
The beetle is so strongly built that I had to use a hammer to drive a pin through 
the specimen I was dissecting. 

Data are as yet lacking to enable us to make a profitable comparison between 
the wings of different insects. On the present occasion I will briefly describe the 
wings of the locust. 


Fig. 1.—Transverse sections of fore wing of locust made at positions shown by 
lines numbered 1 to 6 in Fig. 3. The leading edge of each section is at the end 
near the number. 


The Wings of the Locust. 
A series of sections of the wings of the locust are shown in Figs. 1 and 2. 


It will be noticed that the nervures stand out on both surfaces of the wing 


forming small ridges situated at right angles to the line of flight. The leading 
edge of the whole of the fore wing is very sharp, as is the leading edge of the 
outer half of the hind wing. The strength of the fore wing and the size of its 


nervures are probably partly due to its having to act as a protection to the hind 
wing when the insect is at rest. 


The following illustration, Fig. 3, shows the outline of the Indian locust 
(Acridium peregrinum) with its wings extended. When the insect is at rest, as 
is the case with all grasshoppers, the hind wing is folded up like a fan and lies 
under the front wing. The hind wing of a grasshopper is shown expanded in 
Fig. 4 and folded up in Fig. 5. 


When the wing is partly expanded the folds do not completely disappear 
and the wing membrane forms strong ridges that are partly situated at right 
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angles to the line of flight. It is a question of some interest how far these ridges 
disappear when the wing is fully expanded since I have shown elsewhere that the 
presence of ridges on the underside of the wing is characteristic of the more efh- 
cient soaring birds and other soaring animals (*‘ Animal Flight,’? Chapter XIX.). 
In the case of a specimen of a locust that had been preserved in spirit I found 
that on unfolding the hind wing with my fingers the folds did not completely 
disappear. With an insect that had been kept in captivity till it was moribund 
and then killed the wing could be further flattened. In the case of an insect that 
was killed immediately when caught the unfolding could be carried further till the 
wing was completely flat except at the inner portion near its insertion on to the 
body. But to thus flatten the wing a certain amount of force was necessary and 


* 


Fia. 3.—Outline of locust. The lines numbered 1 to 6 show positions of sections 
of wings shown in Figs. 4 and 5. 


there is room for doubt how far the muscular effort of the insect could produce the 
same result. There can be no doubt that when flapping the outer part of the 
wing is perfectly free of folding. But it seems probable that the inner part retains 
slightly more folding than is indicated in the sections shown in Fig. 2. 


The locust of which I have been treating is known in India as the north-west 
or migratory locust. It is capable of very long flights which all the evidence goes 
to show are entirely by flapping. Gliding appears to be used solely when ths 
insect is descending to the ground. As I learn from an interesting communication 
from Mr. Maxwell-Lefroy, Entomologist to the Government of India, the swarms 
that reach India are ‘‘ probably due to invasions from Arabia, Persia, Afghanistan, 
and Baluchistan, which enter India, breed (in desert portions of Western India), 
and the resulting swarms pass on eastwards to perish when they get beyond the 
sandy deserts in which alone they can exist permanently.’’ A friend of mine 
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informs me that once, when out of sight of land on a P. and O. steamer in the 
Arabian Sea, he passed a swarm of locusts of such extent that it remained im 
sight on two successive days. 

The Bombay locust (Acridium succinctum) has its permanent breeding place 
in forest areas in Western India, where it frequently forms swarms of small 


Fic. 4.—Hind wing of grasshopper (Peecilocera picta) expanded. 


— 


Fig. 5.—Hind wing of grasshopper (Peecilocera picta) folded. 


‘extent. Larger swarms occur at long intervals. In other parts of India it exists 


permanently in forest areas, but only in small numbers which never give rise to 


‘swarms. As a rule the limit of its flight by day is a few miles. It is stated to 


fly for longer distances at night.* 


* See ‘ The Bombay Locust,’’ by H. Maxwell-Lefroy, published in ‘‘ Memoirs of the 
Department of Agriculture in India, Entomological Series,’’ Vol. I., No. 1, April, 1906, printed 


‘by Thacker, Spink and Co., Calcutta. 
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